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EXECUTIVE SUMMARY 
Structure of the Effort 

The advent of a broad scale, low cost remote sensing 
capability Initiated by the successful launch and oper- 
ation of ERTS has given renewed promise and emphasis to 
the development of regional hydrologic planning models 
especially suited to un gaged watersheds. 

This project was instituted to determine the effective- 
ness of a nev. hydrologic planning model .specifically tail- 
ored to remote sensing inputs, and to determine which par- 
ameters Impact most the minimization of errors associated 
with the prediction of peak flow events. The ''capability 
of accurately forecasting peak flow strongly influences 
the sizing and design of civil works and the process of 
general regional watershed planning. The objectives of 
the effort are summarized in Figure 1. 

Peak runoff events in ungaged watersheds were chosen for 
investigation because: 

° The prediction of peak runoff - throughout their 
watersheds - is of significant importance to users 
for the purpose of planning, sizing and designing 


waterworks . 





FIGURE I 


OBJECTIVE OF EFFORT 


DEVELOP WATER RESOURCES 
PLANNING MODEL: 


• Built upon Romote Sensing 
capabilities 

• Suitable for ungoged watersheds 

• Based, for parameters not oc^ssible 
to Rsmote Sansing, upon oommoniy 

available information. 
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® The user community for ungaRed hydrologic planning 
models is significant, including: 

® Federal Agencies involved in water resources. 

® Federal agencies Involved in funding water res- 
ources research. 

® Developed and developing nations. 

® Foreign governmental and international agencies 
responsible for, or oriented towards water res- 
ources development. 

Prom the start, this project espoused the four theses -syn- 
thesized in Figure 2, and described in the following: 

® Of the many phenomena and relationships • underlying 
the behavior of hydrologic processes governing peak 
flow, only a few are dominant for planning models. 
Others can be neglected or, in che extreme, fac- 
tored in as small, constant corrections. 

o The Incluj-.ion of temporal and areal variations of 
watershed characteristics, which have generally 
been averaged or Ignored in hydrologic modeling, 
should result in reduced errors of prediction. 

o The crucial phenomena regulating peak events are, 
to a significant extent, surface-dependent and thus 



FIGURE 2 
THESIS 


• PEAK EVENTS (SUCH AS FLOODS) ARE CAUSED 
BY THE INTERRELATIONSHIP OF A RELATIVELY 
SMALL NUMBER OF WATERSHED PARAMETERS 


• PEAK EVENT MODELS TO BE ACCURATE MUST 
INCORPORATE VARIABLE PARAMETERS 

• REMOTE SENSING TECHNIQUES CAN BE BUILT 
INTO HYDROLOGIC FLOOD MODELS AB INITIO 

• GOOD RESULTS CAN BE OBTAINED WITH 
RELATIVELY SMALL COMPUTER HARDWARE 
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hlghly amenable to remote sensing. For the cases 
where subsurface phenomena play a significant role, 
peak events can still be predicted through a com- 
bination of remotely sensed observables and avail- 
able subsurface parameters. 

® Exploitation of the sensitivities of the physical 
factors underlying hydrologic phenomena can lead 
to significant modeling simplifications from which 
economies in computing hardware will accrue. 

To test the validity of these hypotheses, four investigations 
were undertaken; 

1. A mathematical sensitivity analysis of the hydrologic 
processes cogent to peak watershed outflows was under- 
taken to ascertain which had the greatest influence. 
Further, physical characteristics of basins, such 

as soil type and vegetative cover, were examined 
to determi.ne the sensitivity of peak flow to changes 
in their antecedent conditions. 

2. The construction of a generalized hydrologic plan- 
ning model was initiated along two courses Analog 
computer circuits describing hydrologic processes 
were assembled and applied to examination of the 
sensitivities described above. Also, an analytic 
formulation for peak flow prediction was developed 
and tested against data from actual watersheds. 
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3. A data base of 158 watersheds was collected, con- 
taining physical parameters, rainfall and flow rec- 
ords. The application of this information was two- 
fold - to provide a measure of the validity of the 
models developed, and to permit the delineation of 
areas within the U.S. which are the most amenable 
to modeling from remotely sensed data. Their loc- 
ations are shown in Figure 3. 

4. The role of remote sensing in hydrologic modeling 
was identified and a procedure for its use demon- 
strated. 

The interrelations of these sub-tasks are depicted in Fig- 
ure 4. 


TECHNICAL PROCEDURES 


The Drivers 

Examination of the rates and magnitudes of hydrologic proces- 
ses, indicated that the following processes are important 
to peak flow prediction: 

1) Precipitation 
' 2) Infiltration 

3) Overland Flow 

Secondarily, according to the region under consideration, 
the statistical behavior of antecedent moisture condition 



FIGURE 3 AGRICULTURAL RESEARCH BASINS 




FIGURE 4 PROGRAM FLOW CHART, NAS8-3O530 
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as it relates to peak rainfall events must be Included. 
Rainfall is the principal causative factor defining the 

’ I • ' * 

magnitude of the peak flow. Its important characteris- 
tics are the recurrence statistics, determined by empir- 
ical, correlation of regional rainfall records. 

Infiltration governs the portion of the rainfall which 
contributes to the direct runoff peak. It can be eval- 
uated from watershed soil records, abundantly;- available . 

The overland flow process and channel flow determints 
the timing f the peak. The timing in turn determines 
the rainfall rate and mass for a given recurrence frequen- 
cy, and hence determines the peak flow. The overland flow 
can be modeled from knowledge of the surface characteris- 
tics of the watershed, which are directly amenable to re- 
mote sensing. 

The relative contribution of each of the above factors 
to the peak event varies with the antecedent conditions 
of a specific set of watershed parameters. Rainfall, in- 
filtration, and overland flow are sensitive to these par- 
ameters, and their rates and magnitudes are "driven” by 
the particular mix of the parameters present. Thus, the 
rate, volume, and timing of the peak event are directly 
related to this ret of arivers. ' ■ 
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The phenomena which are the key drivers of peak flow. In 
addition to rainfall statistics, are: 

1. Soli Permeability - high permeabilities mear high 
acceptance of water and smaller runoff mass. 

2. Soil Water Capacity - a soil having a greater water 
capacity will retain more rainfall and produce less 
runoff. 

3. Antecedent Soil Moisture - as soil moisture rises, 
the soil becomes more saturated, slowing Infiltra- 
tion rates, reducing total soil moisture capacity, 
and Increasing the runoff volume. 

4. Slope - flow velocity varies directly but non-lin- 
early with slope. 

5. Surface Friction -- velocity varies inversely and 
non-llnearly with surface frlct:'on. 

6. Drainage Density and Pattern - defines the relative 
distances that water will flow overland and in the 
channel in combination with slope and surface fric- 
tion, and defines concentration time. 

Data acquired from the 158 test watersheds provided the 
basis for partitioning the United States into areas ev- 
idencing either surface-dominated or subsurface-dominated 
hydrologic regimes based on the relative magnitudes of 
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the drivers listed above. An Initial partition of the 
United States into 3 categories of hydrologic regimes 
was made. The three regions are: 

1. Heavily Surface Dominant - Where the percentage 
of rainfall to runoff significantly exceeds the 
percentage of rainfall to infiltration. 

2. Surface Dominant - Where more rainfall runs off 
than infiltrates . 

h 

3. Subsurface Dominant - ^iHierfe more rainfall infil 
trates than runs off. 

It should be noted that the regions which are surface dom- 
inated, and, therefore, most amenable to modeling from 
remotely sensed information are also thos which have his- 
torically experienced the greatest flood damage. The In- 
itial partition is present^ in Figure 5. 

The Analytic Model ^ 

The analytic model was developed modularly. The modules are 

\ 

1. Rain recurrence - by examination of U.S. records, 
a general expression of the form: ^ 

T«2 

i * -i 

(t+d)"^> 







AGRICULTURAL RESEARCH BASINS 100-100,000 ACRES 
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Surface Dominant 
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Vfhere: 

1 ■ rain rate m/sc . 

T « recurrence period, years 
t » duration, hours 

a^,a 2 >oi 2 »d = constants, function of location 

has been developed. The coefficients vary among regions, 
but can be determined for each. By this means, the in- 
tensity of a rainfall of any duration and recurrence, with- 
in any region, can be predicted. 

2. Rainfall Spatial Correction - In the larger water- 
sheds the use of point rainfall in a hydrologic 
model does introduce errors. For accuracy, a 
correction factor must be included of the type: 

Pe * 

Vhere : 

Pg * effective rainfall rate, cm/hr 
P = point rainfall rate, cm/hr 
C * correction factor 

3. Subsurface Abstraction s - based upon analysis of 
various formulations, the Holtan equation for in- 
filtration has been adopted for the time being to 
describe the subsurface precipitation losses. 
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a • GI * (Sa - l)^*^ + if 


Where: 


I * total infiltration rate (or subsurface 
abstraction rate) 

a “ average vep;etative cover factor 
QI * maturity of cover 

Sa = averapje available water capacity = total 

available storage - initial moisture content 

I = cumulative infiltration 

If = final infiltration rate 


4 . 


5 . 


Overland Flow 

A closed form solution for the unit peak flow from 
a simple unitary watershed (Flpjure 6) was formulated. 
This formula is primarily dependent on variables 
which can be obtained by remote sensing techniques. 
The equation takes the form: 


Stex = 





Point-of-FIoodlnp; 

In roost practical applications, the user of a hydrologic 
planning model is Interested not only in the accurate 
value of the peak flow, but in the coordinate where, 
along its length, the stream or channel will actually 
begin to flood. The stream begins to flood when 
the water level equals the height of the banks. 
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I - Average overlandflow 
length, meters 
N» friction coeff. 

S » Slope, meters/meter 
t » Roin duration, hr$. 

T ■ Recurrence Interval, 
years 

= Regional 
' rainfall 
ooefficienrs, 

derivoble 
from rainfall 
records 

g»Ke<,T^* 

L> Channel length, m 

(j)= Channel parameter to 
overland parameter ratio 



eC, r-'t 
(t ♦ df3 




d) (£n)^° 

3600 S^/'° §*/* 



PEAK FLOW FOR SURFACE- DOMINATED 
WATERSHEDS-OVERLAND FLOW CONTRIBUTICH^ 


FIGURE 6 
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Floodin.w does not necessarily have to occur at the 
watershed's outlet; it is a function of channf.l 
shape, slope, and rourhness. These factors are con- 
veniently combined into a single formulation ex- 
pressed in terms of the most easily observable para- 
meter, namely channel width. Ploodinp" bepjins to 
occur when channel width (w ) at any point along the 
channel is less than: 


VJhere: 





% 


» channel length, m 

C = roughness ratio surface/channel 

k = channel geometric correction factor 

Yq = depth of overland flow, m 

* channel width at distance from 
bep^lnning of channel 


Figure 7 supplies an example for typical values of 
the parameters. 


A test of the model's validity was made on a randomly 
selected subset of nine of the 'VRS watersheds. These 
were selected to maxiraize geographic and physical div- 
ersity. A detailed analysis of each was performed to 
quantify those parameters renulred to onerate the mod- 
el. A summary of data collected fo'^ the Coshocton, 
Ohio watershed is presented in Figure 8 as an example. 
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THE EFFECT OF OVERLAND FLOW 
ON FLOOD CONDITIONS 


[Average Depth of Overlond Flow- cm 



, 3/8 c 3/8 ^ 5/8 
Lc b Yq 

)l/4— 

Lq- channel length-m 

§ " roughness ratio 
overlond/chonndl 

k * channel geometric 
correction factor 

Yo* depth of overland 
flow-m 

Wq” channel width 
0 outfolhm 

Example shown : 

k« l/K) §»|.| 


6 9 12 15 

Channel Outfall Width-m 
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RGURE 7 


RGURE 8 

COSHOCTON OHIO ECO-2 


Area » .76 km^ 

Slope ■ .172 m/m 

Shape" Square 

Length of Chofmel - 1491 m 

Drainage Density" 1/510 m/m^ 


*m/m< 


Cover 


Soils 


.039 

(t ^.2) 83 


23% Hardwood Forest 
58% Grossland 
1 1 % Cuitivoted 
8% Miscellaneous 


33% Muskingum silt toam 
19% Keene shallow loom 
17% Keene silt loam 
17% Mixed silt loam 
14% Muskingum stony loom 




The analytic model was applied to each of the watersheds. 


Verification of the Model 

Equally applied were three other hydrologic planning mod- 
els In wide current use - the fetional formula method. 
Cook’s method, and the Soil Conservation Service method. 

The 50-year peak flow statistically determined was used 
to compare the model’s predictions with reality. Figures 
derived are reported in Figure 9. 

The results of thw tests, albeit on this linitod sample, 
shovrcd that a rcnote scnsinf^ oriented planning model 
has the potential for reducing predict ic.i errors that 
occur in conventional models for relatively small un- 
itary ungaged watersheds. The balance of the effort 
centered on identifying the applicability of current 
and near-future remote sensing techniques for developing 
the information required by the model. 

A task was initiated to derive the information required 
to operate planning models from satellite imagery . To 
this end, a 4x enlargement of an ERTS photograph (Figure 
10) of the Chickasha, Oklahoma test watershed vjas examined. 
It was determined that basin area and drainage density 
could be quantified directly from the picture, and that 
surface cover could be identified. The assessment of 
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FIGURE 9 

COMPARISON OF RESULTS 
FOR PEAK OF THE FIFTY 
YEAR EVENT 


1. Dofivme Vt. 

2. Coshocton, Ohio 
Z. Blacksburg. Vo 

4. Oxford Miss. 

6 Fennimore Wise. 

6. CNckasha, Okie 

7 Waco, Texas 
8l Safford, Ariz. 


OgO " m 3/ sec/km 2 

Records 

ECO 

— 

Rattonoi 

SCS 

Cook 

0.95 

0.91 

4.8 

2.14 

5.49 

10.6 

25.5 

17.6 

4.4 

12.6 

1.33 

iOi 

12.7 

7.5 

II. 1 

11.9 

lOB 

7.3 

3.1 

8.4 

ii.d 

12.5 

16.8 

3.5 

13.1 

0.88 

0.08 

3.3 

2.9 

6.44 

13.6 

11.5 

?5.4 

223 

5.7 

6.25 

5.3 

14.4 

15.2 

. 

5.0 

0.87 

0.001 

1.7 

13.7 

3.9 


9. Reynolds, Ohio 
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4x ERTS IMAGE, CHiCKASHA,OKLA. 
R6URE 10 WATERSHED 
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subsurface conditions would require Inference from cover 
data, but this Is beyond the scope of this effort. It 
appears that current experimental techniques for extract- 
ing Information from Imagery and computer- compatible tapes, 
when fully operational, could play a significant role in 
obtaining this information. 

Results and Conclusions 

1. An improved model for the prediction of peak flow 
events has been structured, which is specifically 
designed to take maximum advantage of the data 
and information stream available from remote 
sensing. 

2. The development of the model has been carried to 
the point where the overall framework has been 
constructed and five modules simulating the 
behavior of significant hydrologic processes 
have been developed. 

3. The improved model is considerably more sophls- 
. tlcated than conventional hydrologic planning 

models. In particular, its modules are not 
simply interconnected, but require feedback. In 
spite of this greater complexity, however, the 
model is readily adaptable to analog computation 



with modest amounts of hardware. Preliminary 
sizing shows that the technique can also be 
programmed onto one of the smaller types of 
digital minicomputers . 

The model was exercised — not in its fully 
Interconnected form, but rather in a simplified 
version — to predict the peak runoff from nine 
experimental Agricultural Research Service 
watersheds, selected at random from among a 
set of 158 instrumented and well-described 
watersheds . 

The predictions of the new model in its simpli- 
fied version were tested against: 

a. The predictions from three of the most 
employed contemporary planning models — 
l.e., the rational formula method. Cook's 
method, and the Soil Conservation Service 
method. 

b. The streamgage records of the nine test 
watersheds . 

The results indicate that, v/lthln the range of 
applicability of its simplified version, the 
new model appears to be considerably more ac- 
curate than conventional hydrologic planning 



models. Specifically, In six out of nine of the 
watersheds tested, the new model supplied pre- 
dictions of peak flow for the 50-year event 
fflllnp within error bounds of + 1 ^%, For these 
same six watersheds, conventional models yielded 
discrepancies v;lth respect to the records ranging 
from a minimum of 1.2 to 1 to a maximum of 15 to 
1. For the three reralnlnp watersheds, the 
new model yielded predictions of lesser accuracy, 
the worst being 2 to 1. Reasonable explanations 
for the discrepancy are: 

a. The fact of having oversimplified the model 
by not operating it in its fully intercon- 
nected form. 

b. The three watersheds are considerably more 
complex than the other six, and they need to 
be split into sub water sheds, predicting the 
output from these, then routing all outputs 
through the watershed channels. This 
technique, which appears to be well in hand, 
is proposed for future phases of this effort. 

The appropriate technioues whereby to extract the 
inputs and parameters required by the new model 
from remotely sensed information — whether 
imagery or digital tapes — were explicitly de- 
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flned. Their feasibllitv was identified from 
specific past and onprolnjr ERTS investigator 
efforts. 



